Dynamical evolution of an inverted spin ensemble in a cavity: 
Inhomogeneous broadening as a stabilizing mechanism 
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We study the evolution of an inverted spin ensemble coupled to a cavity. The inversion itself 
presents an inherent instability of the system; however, the inhomogeneous broadening of spin- 
resonance frequencies presents a stabilizing mechanism, and a stability criterion is derived. The 
detailed behavior of mean values and variances of the spin components and of the cavity field is 
accounted for under both stable and unstable conditions. 
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I. INTRODUCTION 

As an extension of traditional cavity-quantum-electro- 
dynamics the resonant coupling of a cavity to an 
ensemble of two-level systems has received considerable 
interest for the past three decades. In particular, by 
the collective effect of N particles, the otherwise weak 
single-particle coupling is enhanced by a factor of >/~N 
0. In experiment, this has allowed for reaching the col- 
lective strong-coupling regime in atomic ionic @, 
and solid-state implementations @-H[ , typically material- 
izing as a normal-mode splitting of the coupled radiation- 
matter system. Ensembles of electronic spins coupled to 
a micro-wave cavity have recently been considered for 
quantum-memory purposes [10h14| |. However, such en- 
sembles usually contain an inherent inhomogeneity of the 
spin transition frequencies, which leads to dephasing of 
the stored information. While spin-refocusing techniques 
reverse this process, it is necessary to understand the dy- 
namical effects and stability of an inverted ensemble cou- 
pled to a cavity-field mode to benefit from such refocusing 
processes. This is the topic of the present manuscript. 
We demonstrate, in particular, that the inhomogeneity 
of the spin transition frequencies is an advantage in the 
sense that it plays a stabilizing role for an inverted ensem- 
ble. Dynamical effects will be examined for both mean 
values and second moments, and a stability criterion for 
an inverted sample is derived. 

Our analysis applies in general for any large collection 
of two- level systems, but for convenience we shall use the 
terminology and notation of ensembles of spin- ~ parti- 
cles. The paper is arranged as follows: In Sec. iHTthe ba- 
sic interaction and decay mechanisms of the spin-cavity 
system is described, and the dynamical evolution of the 
physical system is calculated with emphasis on mean val- 
ues in Sec. IHII and on second moments in Sec. IIVI A few 
experimental diagnostics tools are suggested in Sec. [Vj 
while a general discussion and conclusion of the results 



are given in sections I VII and IVII1 respectively. Some 
mathematical details have been deferred to appendix [X] 



II. EQUATIONS OF MOTION 

We consider an ensemble of N spins coupled to a single- 
mode cavity field, a c , as shown in Fig. [TJ The resonance 
frequency, L)j, of each spin is assumed to be inhomoge- 
neously broadened around a central frequency, (j s , and 
the coupling strength, gj, between individual spins and 
the cavity field may also vary. An external field, ft, 
may be used to drive the cavity field through the left- 
most mirror with field-decay rate, n\ (in the present 
manuscript this driving field is only used for diagnos- 
tics and otherwise left at zero). In the frame rotating at 
the central spin frequency, cj s , the Hamiltonian can be 
expressed as: 
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where A c 

onance frequency u c from w s , and Aj 



is the detuning of the cavity res- 



uij — u> s . The 

Pauli operators with k = — , +, z are used to model 
the j'th. spin. The c-number, j3, represents an external 
coherent-state driving field and is normalized such that 
|/3 1 2 is the incoming number of photons per second. 

Decay mechanisms are taken into account in the 
Markovian approximation of memoryless reservoirs: Cav- 
ity leakage is parametrized by the total field-decay rate 
k = Ki + K2, and the dephasing rate 7j_ = - represents 
the loss of coherence at the level of individual spins with 
a characteristic coherence time r. 



By defining b 
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Sons is given by £i 2 ns = J2j=i 9j, the interaction part of 
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Spin ensemble 



imation (5z ~ 1) in absence of external driving: 




FIG. 1. The physical setup under consideration. A spin en- 
semble is coupled to a cavity field a c , which is subjected to 
decay through the two mirrors with field-decay rates Ki and 
K2- For diagnostics purposes an external driving field /3 may 
be applied giving rise to reflected and transmitted fields, <zr 
and St, as discussed in Sec. IVl 



the Hamiltonian ([1} can be written as h~g cns (a c tf + alb). 
This becomes particularly useful when essentially all 
spins are in the ground state (ai ~ — 1 and ~ 1 

in which case the spin system can be represented by a 
harmonic oscillator - the so-called Holstein-Primakoff ap- 
proximation [Hj]). The resulting formal equivalence be- 
tween quantized fields and collective spin degrees of free- 
dom, and their coupling strength which is collectively en- 
hanced by a factor of yN, have paved the way for u sing 
spin ensembles for quantum information purposes [10l — 
[l4| . The Holstein-Primakoff Hamiltonian is quadratic in 
the oscillator quadrature operators, which implies that 
first and second moments of those operators are described 
by a closed set of equations, also in the presence of inho- 
mogeneous coupling [16| and broadening. The influence 
of inhomogeneous broadening on a spin-cavity system has 
been studied previously under the Holstein-Primakoff ap- 
proximation for ensembles essentially in the ground state 
(l7H2g |. For an inverted ensemble containing 10 two- level 
systems the evolution of the mean values was studied phc- 
nomenologically in Ref. [21j . The present manuscript is 
focused on large inverted ensembles, in which case the 
convenient description of both first and second moments 
under the Holstein-Primakoff approximation is possible. 



III. DYNAMICAL EVOLUTION OF AN 
INVERTED MEDIUM INSIDE A CAVITY: MEAN 
VALUES 

The present section considers the evolution of mean 
values of the cavity field and the spin components for an 
inverted spin state. The calculations assume a resonant 
coupling between the cavity and the spins, A cs = 0, in 
which case the effects under study are strongest. 
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= -( 7 _ L + 2 A J )(a 0) ) +igj(a c ) 
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We note that if (a c ) is real and positive, the second term 

of Eq. (|3|) will drive (a_ ) toward positive imaginary val- 
ues. In turn, the second term of Eq. @ will drive (a c ) 
further along the positive real axis, and the physical sys- 
tem is thus unstable due to the gain provided by the in- 
verted sample. This scenario resembles to a large extent 
a laser, and normal laser operation is initiated when the 
gain medium is able to balance the optical losses of the 
cavity; however, the case under study here differs from 
normal laser operation by the fact that the inverted spin 
medium behaves coherently. Accordingly, a large cavity 
loss (i.e. a large n) is not the only way to counter-act the 
inherent instability, but dephasing due to inhomogeneous 
broadening will also contribute. 

In analogy to threshold conditions for normal laser op- 
eration, a stability criterion can be derived for our spin- 
cavity system by searching for a critical cavity-coupling 
parameter, k c , which allows for a non-zero steady-state 
solution for (a c ) and ((Tj}. Then increasing (decreas- 
ing) solutions versus time are expected when k < K c 
(k > k c ). To this end, consider first Eq. © in steady 

, which inserted into Eq. in 
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steady state leads to: (n + iA cs )(a c ) = V . — 

J "y_L i^^-^j 

We shall restrict ourselves to inhomogeneous broadening 
with Aj distributed symmetrically around zero, in which 
case A cs = is indeed the relevant choice. The above 
equation can be satisfied for a non-zero (a c ) provided 
that k attains the critical value: 
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where we assumed the distributions of gj and A^- to 
be uncorrelated. Furthermore, the continuum limit was 
taken by using the spin-resonance-frequency distribution 
/(A) normalized such that /(A)dA = 1. The char- 
acteristic width r of the inhomogeneous distribution was 
implicitly defined, and the requirement of k > k c for 
stability can be reformulated in terms of the effective co- 
operativity parameter, C: 



q 2 

C=S2BB. < 1. 



(5) 



A. A stability criterion using the effective 
cooperativity parameter 

Consider the following mean value equations, which 
have been derived under the Holstein-Primakoff approx- 



B. Homogeneous broadening 

Even though the main focus of this paper is inho- 
mogeneous broadening, it is convenient to know the ef- 
fects of homogeneous broadening for comparison. From 
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Eq. ((4]) it follows immediately that r = 7j_ in this case 
(/(A) is a S- function). In fact, Eqs. (f2|) and © can be 
reformulated in terms of the effective spin component 
S cS = £f =1 f^\ where f = £f =1 <?f /AT, and the 
inverted spin-state problem is only two-dimensional: 



d_ 

dt 
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-(n + iA cs ) -ig 
igN — 7j_ 
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cfl'V 



(6) 



On resonance, A cs = 0, the eigenvalues of this linear set 
of equations are: 



4(c- i)«r 



(7) 



Clearly, when C < 1 both eigenvalues are negative and 
the inverted spin state with I" — < -"' 1 ' 
solution. 



= (ST) = is a stable 



C. Inhomogeneous broadening 

In order to examine the dynamical evolution of the 
spin-cavity system with analytical methods in the case 
of inhomogeneous broadening, it is convenient to treat 
Eqs. and §3§ in Fourier space. In order to handle also 
exponentially increasing solutions, we re- write the dy- 
namical variables as (a c ) — (a c )e vt and (a^) — (a_)e qt . 

Assume also that (a c ) — (cr^) — when t < 0, which 
indeed presents a mathematical solution to the differen- 
tial equations. Then, at t = we change abruptly the 
cavity-held mean value (a c ) — > a and study the subse- 
quent dynamics. This scenario is governed by a modified 
version of Eqs. © and © taken at resonance, A cs = 0: 



d(a c ) 
dt 

d(d { i ) ) 

dt 
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= aS(t) - (K + r])(a c ) -i^2g 3 (a 



(8) 



- (7x + V + (° { - > + *9j (Sc) • (9) 



The latter of these can be integrated formally: (a^P (t)) — 

igj Jo e~ i ' 1±+v+lA ^ t ~ t ' ) (& c (t'))dt' , which in turn can be 
inserted into Eq. ©: 

" ! " ! - aS ( t )-( K + r,)(a c ) + f K(t-t')(a c (t')), (10) 



dt jo 

where K(t) = J2f=i 9j e~^ +7J+lA ^ t . Now, by defining 

the positive-time version of K by K + (t) — K(t)-6(t), 
where 8(t) is the Heaviside step function, and by remem- 
bering that (a c (t r )) — when t' < 0, the above inte- 
gration can be extended to plus/minus infinity. Using 
the Fourier transform, (a c (cj)) = j_ (a c (t))e tut dt and 

(5c(*)> = ^ SZo^^e-^dLJ, we find: 



(a c (w)) 
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(11) 



10 3 



10-' 



1^ a ° D -H D a°5.nnnn qan 



10 1 



10" 



qQO °- 20 ^2PQ PQQOi 




FIG. 2. (Color online) Cavity-field decay versus time when 
the inhomogeneous broadening is Lorentzian [dashed lines, 
from Eq. (|15|) ] and Gaussian [symbols, numerical simula- 
tion]. In all cases, g cns = 2r and 7x = 0, while k is varied 
such that C attains the values of 0.05 (magenta tip-up trian- 
gles), 0.2 (red circles), 0.5 (green squares), 1 (blue diamonds), 
and 2 (black tip-down triangles). The vertical axis contains 
(Xc) = ((So) + (&t))/V2. In the range Ft > 2.5 the slope 
of the Gaussian-broadening decay curves are compared to a 
numerically determined value (solid lines) as discussed in the 
text prior to Eq. (|17[1 . The solid curve through the magenta 
tip-up triangles is given by Eq. (|18[1 . 



The Fourier transform K + (ui) can be expressed in the 
continuum limit as: 
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A — to — i(j± + 77) 



(12) 



1. Lorentzian broadening 

For a Lorentzian broadened spin ensemble with / (A) = 

A*+w*/v where w is thc FWHM (full width at half 
maximum), the characteristic width of Eq. (|4]) becomes: 
r = y +7-L- Furthermore, Eq. (fT^j) can be written (using 
the residue theorem): 



9cns 



T + Tj 



and inserting this result into Eq. (1TTT) leads to: 
a(iu! — r — rj) 



(a c (u)) = 



(u> - i[X+ - r)])(u - i[X- - t]})' 



(13) 



(14) 



where A± are the solutions given in Eq. (J7J) . The inverse 
Fourier transform is now invoked, leading to (t > 0): 



<6c(*)> 



(A + +r)e A + t -(A_+r)e 



A-i 



A- 



(15) 



and (a c (t)) = when t < 0. This expression is inde- 
pendent of 77 as it should be; however, for the Fourier 
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transform (a c (w)) to exist, the condition rj > A + must 
be fulfilled, which in fact also ensures that both poles in 
Eq. (fT4|) reside in the lower complex half-plane. 



2. Gaussian broadening 

For a Gaussian broadened spin ensemble with /(A) = 
-pi — e _A / 2<t a where cta is the standard deviation 

v2tt(7a 

of the distribution (connected to the FWHM by w — 
a a \/8 ln(2)) the characteristic width ((3]) is given by: 

r = ^/f ^f^j > where the complex error function is given 

by w(z) — e~ z erfc(— iz) with erfc( • ) being the comple- 
mentary error function [23 and z = % 2 1 - . Equation (fT2l 
reads in this case: 

£+( w ) = J*iL w{S)t (16) 

V 2 CTA 

where z = " + ^/^~ +7 ^ ■ If i s n °t possible to write a gen- 
eral analytic expression for the inverse transform (a c (t)}; 
however, we shall calculate (a c (t)) by numerical integra- 
tion of Eq. (|A9p with a real, non-zero X c — \/2a as 
initial condition at t = 0, and limiting cases will be 
compared to analytical estimates. Such numerical sim- 
ulations are shown (with symbols) in Fig. [2] for various 
values of the effective cooperativity parameter C, and 
comparison to the Lorentzian-broadened case is made 
(by dashed lines, maintaining k and V). The following 
points can be noted: (I) the initial decay seems simi- 
lar for Lorentzian and Gaussian broadening, (II) in the 
long-time limit for Gaussian broadening the decay seems 
to be exponential, and (III) when the coupling is weak 
(C -C 1) the curves for Gaussian broadening appear to 
have a significant quadratic content when plotted on the 
logarithmic vertical scale. 

The single-exponential parts of the decay curves (with 
rate A) correspond to the poles of Eq. (Tnj) . i.e. solu- 

2 

tions {uj — iX) to the equation k — iuj = y^^^w{z) 
taking rj = 0. Provided that |A| S> oa for the initial 
fast decay, we take advantage of the series expansion, 
w(z) sa -j=^ when \z\ >• 1, and reach the condition: 

(A + k)(A + 7_l) = <7g ns . This is exactly the eigen- value 
equation for the homogeneous system of equations ([5]), 
and the solution A is equal to A_ in Eq. ([7]) with r = 7j_ ■ 
From a physical perspective, the narrow feature of the 
Gaussian broadening cannot be resolved on the initial 
fast time scales. The long-time limit of the decay for 
Gaussian broadening is compared in Fig. [5] by solid lines 
to the rate A found by locating numerically another pole 
of Eq. (fTTj) . In the vicinity of the stability threshold, 
k k c such that |A| <C cta, the value of A can be approx- 
imated by using the series expansion, w(z) ss 1 + =j£ — z 2 



when \z\ <C 1, leading to: 

Finally, the weak-coupling limit, C <C 1, can be calcu- 
lated directly from Eqs. flSJ) and ©, provided that k is 
faster than the remaining dynamical processes. In a first 
approximation, (a c (t)) — ae~ Kt , since the spins will con- 
tribute little due to the low coupling. Secondly, during 
the initial decay of the cavity field, each spin component 
acquires a small value: — , which is derived by 

integrating the second term of Eq. ([9]) ; the first term can 
be neglected on this fast time scale. Thirdly, after the 
initial cavity decay, the spins evolve freely due to the low 
coupling: Jj(cr_ ) = — (7_l + iAj)(<r^ and the cavity 
field follows the spins adiabatically in this regime: 

■ N 2 

<a c (t)) » -£5X*«> = ^e-i°* t2 -^, (18) 

where the continuum limit of the inhomogeneous fre- 
quency distribution was taken in the last step. Alter- 
natively, when C < 1, Eq. (|11[) can be approximated: 
(aJuj)) » -^ + 7 ^Ls./T^l ~ + /|HM 

\ °V // K — IUJ (K — tUJ)^ V 2 (TA K — IU K V 2 CTA 

where the second step considers only the low-frequency 
parts of the second term (w{z) varies on the frequency 
scale of a a Ck). The inverse Fourier transform of this 
approximated (a c (oj)) is ae~ K * plus the term found in 
Eq. (fl~8|) . The lower curve (magenta tip-up triangles, 
C = 0.05) in Fig. H] follows E q- © to a large extent. 

IV. DYNAMICAL EVOLUTION OF AN 
INVERTED MEDIUM INSIDE A CAVITY: 
QUADRATIC MOMENTS 

The calculation of the dynamical evolution of mean 
values in the preceding section presents one of the main 
results of the present manuscript. However, we wish 
to back up these mean-field results by a calculation of 
second moments — a mean-value stabilized spin-cavity 
system would be of less relevance if e.g. the variance of 
the spin components and the cavity field increased with- 
out limits. Such an unlimited increase will also render 
the spin-cavity system inapplicable for quantum-memory 
purposes. 

The case of homogeneous broadening is treated analyt- 
ically while inhomogeneous broadening requires numer- 
ical treatment. In any case, the calculations follow the 
general procedure outlined in appendix [S] 

A. Homogeneous broadening 

Assume that Aj = and gj — g for all spins. For 
an inverted spin sample on resonance with mean values 
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S X Sy 0, S Z 



N , and A r 



troduce the vector of second moments, x = [(8X£) (6P ( 



0, we in- 

2\ /ip2\ 



(SS 2 ) (SS 2 ) (5S X SP C ) (SS y SX c )} T . Following Eq. (XI]} 
in the appendix, they obey the following set of coupled 
equations, ^ = Qx + r, where 
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, r = 


2 7± iV 





-V2gN 
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V2 





-(/s+7x) 










V2gN 








9 

V2 





-(k + 7_l)_ 








(19) 



In fact, an inhomogeneous distribution of the coupling 
constants, gj, can be incorporated in the above equations 
by merely replacing S x — > S x s , S y — > Sy S , and 5 — > g. 
The matrix Q has three doubly-degenerate eigenvalues. 
Two of these are given by A = 2A±, i.e. by twice the 
values found in Eq. J7]), and the third one is A = — (k + 
7_l). Hence, the same condition C < 1, ensures that both 
the first and second moments are stable and converge to 
their steady-state values. Solving ^ = Qx+r = 0, these 
read: 

(5X 2 C ) = (SP 2 ) = \ 1 -^^, 



(SSf 2 ) = (SSf 2 ) =N- i+ ° K+r 



(SSfSP c ) = (SS c v n 8X c 



1-C ' 

2g c 



(20) 



2 (K + r)(i-c) : 



where r = 7 j_ for homogeneous broadening. We note 
that the levels of \ and N correspond to the variance of 
the minimum-uncertainty states for the cavity field and 
the collective spin, respectively. When approaching the 
stability point C — > 1 from below, the variances diverge. 



B. Inhomogeneous broadening 

For the case of inhomogeneous broadening we use nu- 
merical simulation of Eqs. (|A9|) and (|A11|) for calculating 
the dynamical evolution. We use only a single value for 
g m but choose either a Lorentzian or Gaussian shaped 
distribution of A m . As a starting point, all spins are 
prepared in the inverted coherent state being slightly 



displaced: S. 



(to) 



(m) 



0, and Si m) = 9Si m) 



where = 10 -3 , and the cavity is prepared in the vacuum 
state. Leaving the spin-cavity system to evolve from this 
initial state, we study as function of time a representa- 
tive set of mean values and variances: S x — J2m=i 

Pc (6S%) = EtL=i(^i m) SS x n) }, and (<5P C 2 ), the results 
have been plotted in Fig. [3] 

Panels (a) and (b) of this figure show how the vari- 
ances, (5S 2 .) and (SP 2 ), increase from their initial val- 



ues of N and i, respectively. As can be seen, in the 
stable region with C < 1 these variances converge to a 
steady-state value while for C > 1 the curves increase 
without limits. The solid lines correspond to a Gaus- 
sian distribution while the dashed lines correspond to a 
homogeneously broadened sample with r = 7 j_, which 
coincides with the simulations for a Lorentzian broad- 
ened sample with 7 j_ = and V = ^. At the same 
time, the mean values of S x and P c have been plotted in 
panel (c) and (d), which confirm that solutions increase 
or decrease versus time when C > 1 or C < 1, respec- 
tively. We note that the features and interpretation of 
these graphs are very similar to those of Fig. [2] only the 
initial state is different in the two figures. In order to 
show the time scale of the dynamical evolution of vari- 
ances, the deviation of (SS 2 (t)) from its asymptotic value 
of (SS 2 (00)) has been shown in panel (e) relative to the 
entire dynamical range, i.e. the vertical scale is the ra- 

t-: R(t) = ;g:j;:gfg; - Noting that in panel (e) 
the horizontal axis spans only half the time as compared 
to panel (c), it can be seen that the variance (SS 2 ) ap- 
proaches its asymptotic value approximately twice as fast 
as the decay of the mean value S x toward zero. This is no 
surprise for the homogeneous or Lorentzian case since we 
already observed that the three characteristic eigenvalues 
of the problem, 2A+, 2A_, and A+ + A_, relate closely to 
the eigen values of the mean value equation • The sim- 
ilarity of the solid lines in panels (c) and (e) demonstrate 
that this holds qualitatively also for the case of Gaussian 
broadening. Finally, it can be observed from panels (a) 
and (b) that in the case of Gaussian broadening (solid 
lines), the variances (6S X ) and (SP 2 ) converge to val- 
ues which are slightly higher than those give by Eq. (|20)) 
when simply inserting the corresponding values for k, T, 
and C (dashed lines). The ratio of solid-to-dashed lines 
in panels (a) and (b) have been shown in panel (f) with 
circles and diamonds, respectively, and varying values for 
the ratio <? e ns/r have been examined. We conclude that 
Eq. ([20|) is not accurate for a Gaussian inhomogeneous 
distribution although the qualitative features remain. 
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FIG. 3. (Color online) Panels (a-e) show various dynamical 
parameters versus time with fixed g cns — 2F and varying k, 
such that C attains the values 0.05 (magenta), 0.2 (red), 0.5 
(green), 1 (blue), and 2 (black), represented in the order from 
the the lower to the upper sets of dashed and solid curves 
in all plots. Solid lines have been calculated for Gaussian 
inhomogeneous broadening (■y± — 0), while the dashed lines 
correspond to both Lorentzian broadening (7j_ = 0) and ho- 
mogeneous broadening (-y± — V). (a) Excess variance of S x 
relative to the coherent-state value of N. (b) Excess vari- 
ance of P c relative to the coherent-state value of i. (c) The 
mean value of S x normalized to S x (t = 0). (d) The mean 
value of P c in arbitrary units, (e) The relative deviation 
of (5S%) from its steady-state value, R(t), see the text for 
details, (f) The Gaussian-broadened steady-state excess vari- 
ance [{6§2(oo)) /N-l] (circles) or [2(<5P c 2 (oo))-l] (diamonds), 
relative to the corresponding values from Eq. (|20[) . The results 
are very similar for the different coupling strengths g cns = 3r 
(red), g cns = 4P (green), and g cns = 5P (blue). 



EXTERNAL PROBING OF THE SPIN 
SAMPLE 



non-inverted sample and also allows for assessing the ef- 
ficiency of the spin-inversion process. Assuming the cav- 
ity to be resonant with the spins, A cs = 0, the follow- 
ing mean- value equations are valid under the Holstein- 
Primakoff approximation: 



N 



at 
at 



= - K (a c ) -i^g^) + V2^P (21) 

3=1 



-(7 ± +iA i )(o-i J ' ) )+ip ft (a c ), 



(22) 



where p = 1 for an inverted sample and p = — 1 for a non- 
inverted sample. By applying a monochromatic external 
field, /3(t) = /3oe~ aAo *, the cavity-field mean value can be 
shown to be: 



<6c(t)> 



P9cn S 



r. 



o /(A)dA 
oo 7_i_+i(A— A B ) 



(23) 



In fact, this is a particular solution to the differential 
equation and we assume that the homogeneous solution 
has relaxed to zero; this relaxation process was the topic 
of Sec. IIII1 and for an inverted sample (p = 1) the calcu- 
lations only make sense if the stability criterion is met, 
C < 1. The integral in the denominator of the above 
equation is equal to [r — iA c ] _1 for Lorentzian broaden- 
ing and equal to y^ w ^ with z c — A ^^ y± f° r Gaussian 
broadening. When the driving is resonant, A c = 0, the 
integral is equal to for any (symmetric) distribution 
according to Eq. 

Now, the reflected and transmitted fields relate to 
the cavity field by [23|: (Sr) = y/2Ki(a c ) — (3 and 
(ax) = v / 2k2(q, c ). This enables a calculation of the com- 
plex reflection and transmission coefficients, r = and 

t = , respectively. Selected examples have been plot- 
ted in Fig. U for a Lorentzian inhomogeneous broadening 
(a Gaussian profile presents qualitatively similar results). 
When C increases beyond unity, the bare-cavity trans- 
mission spectrum is significantly modified by the pres- 
ence of the spin ensemble and the well-known normal- 
mode splitting occurs Note that C = 1 corresponds 
to the case where the transmission coefficient is reduced 
from unity to one half (for a symmetric cavity). We 
also note that the transmission spectrum exists for an 
inverted sample when C < 1 (the dotted curve exempli- 
fies this) and that the transmission coefficient may exceed 
unity due to the inherent gain of the inverted sample. 

A particular relation, which is useful for a simple esti- 
mation of the effective cooperativity parameter, is given 
by the connection of C to the values of r and t for any 
(symmetric) distribution on resonance (A e = 0): 



The linear response of the spin ensemble can be probed 
by applying a weak, external field j3 and measuring the 
reflected or transmitted field as depicted in Fig. [T] Such 
a measurement enables the determination of C from a 



pC = 



Kl — K2 



= 1-^.1. (24) 

1 Kj t 



The fact that the reflection and transmission coeffi- 
cients may exceed unity clearly demonstrates that the 
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FIG. 4. (Color online) The real part of the transmission co- 
efficient versus external driving frequency. The calculations 
assume a Lorentzian profile, Ki = «2, k = 10F, and g cns is 
varied in order to obtain the C-values marked on the curves. 
Solid lines are non- inverted (p = —1) while the dotted line for 
C = 0.2 corresponds to an inverted sample (p = 1). 

excitation level of the spin ensemble must account for the 
energy balance. This fact is disguised by the Holstein- 
Primakoff approximation, but it is possible to estimate 
the effect in a mean-field theory on resonance (A cs = 0, 
see the appendix for details): 

9S Z _ 4 P ffg ns |q c | 2 

~~dt ~ f • (25) 

We stress that this holds also for a non-inverted sam- 
ple (p = —1), in which case energy quanta leak from 
the cavity into the continuous spin ensemble with a 

2 

rate proportional to a squared matrix element, ~ 
times the density of states, ~ resembling the usual 
Fermi's-Golden-Rule expression for decay of a quantum 
system due to the coupling to a broad-bandwidth reser- 
voir. In order that this de-polarizing effect is kept small, 
the duration T of the external driving must be short 
enough so that l^r)? 1 -C N. This is equivalent to: 

n ph < £ ' (1 8|pC| 2 ' N > where n ph = \P\ 2T is the total 
number of photons supplied by the external driving field 
during the experiment. Clearly, for an inverted sample 
approaching the point of instability, pC — > 1, the allowed 
number of photons decreases significantly below N. 



VI. DISCUSSION 

The stability criterion of Eq. ([5]) and the dynami- 
cal evolution of the spin-cavity system below and above 
the point of stability present the main result of this 
manuscript. The results of Sec. IIV| in particular panels 
(a-d) of Fig. [3l demonstrate that the stability criterion 
refers to both the mean values and the second moments. 



This follows naturally from the fact that the same ma- 
trices govern the linear sets of equations for the first and 
second moments, as shown in the appendix. 

Understanding the free evolution of an inverted spin 
ensemble in a cavity is of high importance for spin- 
refocusing techniques. Such refocusing could improve 
spin-based quantum memory protocols in cavities. How- 
ever, the storage and retrieval part of such a protocol [l3[ 
would typically be implemented in the strong-coupling 
regime, g cns 3> K, T, i.e. with C 3> 1, and the ability 
to tune the value of C during the experimental protocol 
would then be necessary. We also note that the spins 
can effectively be decoupled from the cavity field by a 
large detuning A cs . The discussion after Eq. ([T5jl can be 

stated more generally as (a c (u>)) « K+i{ £ a _ u) + 
when |K + iA cs | is large compared to the frequency width 
of K + (uj). This equation reflects the fact that the cav- 
ity field follows adiabatically the evolution of the (ef- 
fectively) uncoupled spin ensemble (the second term de- 
pends on oj through K + (uj) only), which in turn is largely 
given by the Fourier components of /(A) through the re- 
lation (fT2|) . Note that a broad and smooth distribution 
/(A) is required in general, coupled or uncoupled, if a 
fast relaxation of both the cavity field and the spin com- 
ponents is desired. 

The diagnostics tools presented in Sec. [V] have been 
derived for a perfectly polarized spin ensemble (p = ±1). 
However, as exemplified in the appendix by using a suit- 
able sub-ensemble distribution, we may argue that the 
results of Sec. [V] hold for a non-perfect polarization also, 
— 1 < P < 1. The Holstein-Primakoff approximation cor- 
responds to keeping the collective spin vector within the 
linear region around the north or south pole of a collective 
Bloch sphere. Relaxing the need for perfect polarization 
corresponds to merely reducing the radius of the collec- 
tive Bloch vector. We note that the important equations 
include p and C in the combination pC , and since C cx N 
it is reasonable that a non-perfect spin polarization is ac- 
counted for by this combination. This argument holds 
also for the stability criterion of Eq. 



VII. CONCLUSION 

We have demonstrated that inhomogeneous broaden- 
ing is a stabilizing mechanism for an inverted spin ensem- 
ble coupled to a cavity. A stability criterion was stated 
in Eq. ([5]) , and if this criterion is met the transverse spin- 
component mean values relax toward zero while the vari- 
ances of these spin components reach finite values. This 
holds simultaneously for the mean values and variances 
of the cavity field. 

The details of the spin-cavity dynamics was discussed 
for a Lorentzian and a Gaussian inhomogeneity in the 
spin-resonance frequencies. In particular, the time scale 
of the relaxation process is well understood, and fast re- 
laxation requires a broad and smooth inhomogeneity. 
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Appendix A: Equations for first and second 
moments using a sub-ensemble discretization 



The present appendix describes how the Holstein- 
Primakoff approximation is applied to establish numeri- 
cally tractable, linear equations of motion for the mean 
values and second moments of field and collective spin 
quadrature operators. The system Hilbert space is infi- 
nite dimensional even for a single oscillator mode; how- 
ever, the number of equations we need to solve scales only 
linearly and quadratically with the number of modes for 
mean values and second moments, respectively. In or- 
der to handle inhomogeneities numerically, the ensem- 
ble is divided into M sub-ensembles, M.\,M.2, ■ ■ ■ ,Mm, 
which can each be regarded as homogeneous with cou- 
pling strength g m , spin resonance frequency A m , and 
containing N m spins for m = 1, . . . , M. We assume that 
all spins reside in the ground or in the excited state, such 
dj) 



J2m " z ~ i^m, and that the dynamical 



that 

variables of the spin-cavity system are described by the 
operators: 



X c = 



Pr 



a c + at 



V2 ' 
-i(a c - at) 



V2 ' 



(Al) 



M,-, 



3f> = -i5>S»-*W) 



The X c and P c operators describe the quadratures of the 
cavity field with [X c , P c ] = i, while the components 
correspond to twice the total spin in each sub-ensemble 



(m) 



with [Sp,SZ n) ]=2ie j kiS l 

The strong spin polarization ensures the constant com- 



mutator: [Sx" 1 ^ S; 



V 



2iS. 



(m) 



±2i/V m , and validates 



the ensuing simplified Heisenberg equations of motion in 



the Holstein-Primakoff approximation 
dX c 



dt 



dP c 
dt 

dsi m) 

dt 



K X c + A cs P c -Y,^S^+F Xc 

m * 

K P c ~A cs X c -Y,^Si^+F Pc7 

rn V 2 

7 ^ m) - A m S^ - s/2g m S^P c 



(A2) 
(A3) 



f c 



(A4) 



dS„ 



(m) 



dt 



-1±S^ + A m S^ - V2g m SWx c + F s 



(A5) 



Note, the external driving is assumed to be absent, ,0 = 0. 
The last term in each equation is a Langevin noise op- 
erator, the properties of which follow from the quan- 
tum Langevin equations of a damped harmonic oscil- 
lator 24j. For instance, the preservation of commuta- 
tors require that [F Xc (i), F Pc (£')] = 2n-i5{t-t') and 
[F sS n) (t), F s (m) (*')] = 4rf±N m ■ iS{t - t'). 

Arranging the field and spin operators in a column vec- 
tor y with 2M + 2 components, we can write the coupled 
Heisenberg equations of motion in the compact form: 



|y = My + F, 
where the driving matrix M is given by: 



(A6) 



M = 



A B« B( 2 ) 



CM DM 
CM 

C (' M ) 





D M 






d( m ) 



(A7) 



with 
A 



— K 

-A cs 



9m 
V2 



9m 





Q(m) _ 
£)( m ) = 



-V2g m 3; 



(m) 



-V2g m si m) 




(A8) 



"7-L 



-A m 

-7-1 



Inserting y = y + Sy, where (Sy) = 0, into Eq. 
yields the mean value equation for y: 

dy 



(A9) 



This is the equation solved in our numerical mean field 
analysis, and it is the eigenvalues of the matrix M which 
govern the stability of the solutions. The Heisenberg 
equations for the fluctuations around the mean values 



d_ 

Of 



Sy = M_5y + F, 



(A10) 
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are operator valued, and to investigate the fluctuations 
numerically we introduce the (2 + 2M) x (2 + 2M) 
covariance matrix 7 = 2Re{(5y • <5y T )} with elements 
Iki = C(y k ,Vi) = 2Re{(Sy k Syi}}, which represent the 
quantum correlations between any two of the 2 + 2M rel- 
evant operators, yu and yi. In particular j kk = 2Var(j/fc). 
The special form of 7 leads to its time derivative: 



^ = M 7 + 7 M T + N, 

at 



(All) 



where N is related to the correlation of the Langevin 
operators by NS(t - t') = 2Re{(F(i)F(i') T )}, which for 
reservoirs at zero temperature amounts to: 



N = 



V 
U« 





U( m ) 



... 










Aj ± N„ 



V = 



2k 
2k 



We observe that the eigenvalue spectrum of the matrix M 
also accounts for the stability properties of the covariance 
matrix, and thus the second moments of collective spin 
variables. 

The sub-ensemble grouping of spins serves two pur- 
poses. Most importantly, it enables the application of 
the Holstein-Primakoff approximation, which results in 
linear coupled equations for the first and second moments 
of effective oscillator quadrature operators. This signif- 
icantly reduces the number of dynamical variables ac- 
counting for the full quantum state to (2M + 2) mean 
values and (2M + 2) 2 second moments with M being the 
number sub-ensembles. For the validity and accuracy of 
our approach one should ensure that M is large enough 
to adequately represent the inhomogeneous broadening 
of the spin ensemble, i.e., the frequency spacing must be 
sufficiently small to avoid discretization errors such as 
artificial revivals of the spin state, while still treating a 
sufficiently large number of spins to render the Holstein- 
Primakoff oscillator description valid. 

Although we treat the collective S z m ^ operators as 
constants equal to ±N m , the validity of the Holstein- 
Primakoff approximations merely relies on their mean 
values being much larger than their quantum fluctua- 
tions. Our analysis will thus also apply for partly po- 
larized samples, and we may revisit the full Heisenberg 



equations of motion in order to determine if they change 
in time due to the coupling to the cavity field. To this 



end, consider the time derivative of £>; 



(m) 



(S z m) ): 



dS : 



(m) 



dt 



(A13) 



In the case of a seeded cavity, studied in Sec. E] the spin 
and field operators have finite mean values, and the above 
expectation values approximately factor: (<Sj!™^a c ) ~ 
S+ ^a c , etc. From the mean value equation 



dS 



(m) 



= -( 7X + iA rn )S im) + ig m Si m ^a c , (A14) 



dt 



(m) 



(A12) 

we can adiabatically eliminate the spin variable: S_ 
%9r " a ° c t° a good approximation when a c varies slowly 



7A+»A„ 

(A cs = 0). Setting S z — pN m in this expression, and 
inserting the result into Eq. (|A13I) , the change of S z be- 
comes: 



dS z s—^ dSi ^ 1 | 2 9mNm 



dt 



dt 

m 

r° /(A)dA 



(A15) 



The second equality is valid for a symmetric sub- 
ensemble distribution, while the third equality assumes 
the continuum limit of the sub-ensemble description of 
the actual inhomogeneous distribution. Using Eq. §4§ 
leads to the Fermi's-Golden- Rule-like expression (|25[) . 
which both explains the field loss due to absorption by 
the non-inverted spin ensemble and the gain obtained 
due to stimulated emission by the inverted sample. We 
note that in the absence of coherent driving, mean val- 
ues of the cavity field and the spin components vanish 
in steady state, and the product term in Eq. (IA13[) is a 
combination of the second moments obtained by solving 
Eq. (jAllj) . For the special case of an inverted homo- 
geneous or Lorentzian spin ensemble, these second mo- 
ments are given by Eq. (|2"U)) . and the rate of change in S z 

when p — 1 and zero when 



becomes ^f- — 
P = -l. 



(«+D(i-c) 
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